Spectral fluctuations and 1/f noise in the order-chaos transition regime 
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Level fluctuations in quantum system have been used to characterize quantum chaos using ran- 
dom matrix models. Recently time series methods were used to relate level fluctuations to the 
classical dynamics in the regular and chaotic limit. In this we show that the spectrum of the system 
undergoing order to chaos transition displays a characteristic /~ 7 noise and 7 is correlated with the 
classical chaos in the system. We demonstrate this using a smooth potential and a time-dependent 
system modeled by Gaussian and circular ensembles respectively of random matrix theory. We show 
the effect of short periodic orbits on these fluctuation measures. 
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Quantum chaos, the study of quantum analogues of 
classically chaotic systems, is characterized by the fluc- 
tuation properties of the spectrum of its Hamiltonian op- 
erator. For quantum systems with regular classical dy- 
namics the spectral fluctuations are Poisson distributed 
0, i.e, the eigenvalues tend to cluster together. On the 
other hand, one of the remarkable results established by 
Bohigas et al is that the level fluctuation properties of 
quantum systems, whose classical limit is chaotic, are 
identical to those of an appropriate ensemble from ran- 
dom matrix theory (RMT) [2J . This is the level repulsion 
regime where the eigenvalues tend to repel one another. 
In this sense, characterizing quantum chaos in terms of 
presence or absence of level repulsion requires invoking 
the spectral properties of random matrix ensembles. Re- 
cently in analogy with time series, a method has been 
proposed to characterize spectral fluctuations using in- 
herent properties of the spectrum 0. If the eigenvalues 
of Hamiltonian operators could be thought of as a time 
series and its index the time in some units, then the meth- 
ods of traditional time series analysis can be applied to 
it. It was shown that the ensemble averaged power spec- 
trum (S(f)) of the fluctuations in the cumulative level 
density, goes as 1/f or 1/f 2 depending on whether the 
system is classically chaotic or regular |4| . This work also 
showed examples of atomic level sequences displaying 1/f 
noise. Thus, atomic levels join the host of other systems 
and phenomena that display 1/f noise lending strength 
to the well known cliche that 1/f noise is ubiquitous in 
nature [||. 

In this paper, we study the transition from regular- 
ity to chaos in the mixed systems. In such systems the 
regular and chaotic motion coexist and this is a generic 
feature. For instance, the entire class of atoms in strong 
fields and the range of problems involving atoms in time- 
varying fields belong to this class. We study a smooth 
Hamiltonian system, the quartic oscillator and a time 
dependent system, the kicked top. In both these sys- 
tems, a single parameter that controls the classical chaos 
can be varied to get a smooth transition from regular to 
predominantly chaotic dynamics. It is well known that 
the level fluctuations in these systems can be modelled 



by RMT @. From an RMT point of view, these mod- 
els possess symmetries (time-reversal invariance without 
spin-1/2 interactions) such that the kicked top is part of 
the circular orthogonal ensemble whereas the coupled os- 
cillator falls in the Gaussian orthogonal ensemble (GOE) 
of RMT [6j . Based on the numerical evidence from these 
models we show that (£(/)) cc /~ 7 , where 7 depends on 
the degree of their classical chaos. This correlation be- 
tween 7 and the classical chaos parameter is established 
using semi-empirical level spacing distributions studied 
in the context of RMT to model the transition region. 

In semiclassical systems of the type we consider here, 
periodic orbits via the Gutzwiller formalism play an im- 
portant role in determining the quantum spectrum [fj. 
As pointed out by Berry J] , the properties of the spec- 
trum on a scale of mean level spacing are determined by 
long time period orbits and in a sense both of them dis- 
play universality; the RMT type universality in the spec- 
trum and classical universality embodied in the Hannay- 
Ozorio sum rule 8]. However, long range spectral prop- 
erties are determined by the short time periodic orbits 
which are system specific and are not universal. This 
manifests itself in the power spectrum as deviations from 
the / 1 behavior. We show the effect of short time peri- 
odic orbits in the coupled oscillator, where scarring or the 
density enhancements in the vicinity of certain periodic 
orbits is a prominent feature due to these orbits. 

The Hamiltonian for the coupled quartic oscillator is, 



Hi = p\ + p\ + q\ + q\ + a q\q\. 



(1) 



This system is classically integrablc for a = 0, 2, 6 and the 
phase space is predominantly chaotic for a > 6. This has 
been extensively studied as a model for chaotic dynamics 
in a smooth potential |lfjl ITlj . The Hamiltonian is quan- 
tized by solving the corresponding Schrodinger equation 
in the basis of the eigenfunctions corresponding to a = 0. 
Then, the matrix elements of Hamiltonian operator Hi 
is computed in this basis and the matrix of order 13000 
is diagonalized to obtain about 2000 converged eigenval- 
ues. This system possesses C^ v point group symmetry. 
We symmetry decompose the spectrum and study the 
levels from only the Ai irreducible representation. 
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The quantum top 01 * s characterized by an angu- 
lar momentum vector J, whose components (J x , J y , J z ) 
obey the usual commutation relations and J 2 = j(j + 

l),j = i,l, |, , is conserved. The dynamics of the 

top is governed by the Hamiltonian |l2|. H^if) = § J y + 
j]Jz J2 n S(t — n). The first term describes the precession 
around the y axis with angular frequency 7r/2 and the 
second term kicks in periodically with (5-function kicks of 
strength k. Each kick can be thought of as an impul- 
sive rotation about z-axis by an angle kJ z /2j. The time 
evolution operator in between consecutive kicks is, 



r = t m> ( -i^jjf) ex P (~*| J a 



(2) 



If are the eigenvalues of U, then the power spectrum 
is computed from the quasienergy (p. In the limit j — > oo, 
we can derive a classical map whose dynamics depends 
on the parameter k |l2|. At k = it is integrable and 
becomes increasingly chaotic for k > 0. 

We denote the eigenvalues of the appropriate opera- 
tor described above by Ei, i = 1,2, . . . ,n + 1. The in- 
tegrated level density, that counts the number of levels 
below a given E, can be decomposed into an average 
and an oscillating part, N(E) = ~N(E) + N osc (E). In 
order to compare the fluctuations from various systems 
it is customary to unfold the spectrum by a transforma- 
tion Xi = N(Ei), such that the mean level density of the 
transformed levels is unity. All further analysis is carried 
out using the sequence {A}. For instance, the spacing is 
Si = Aj+i — A-j, i = 1, 2, . . . , n. In this paper, we will work 
with the statistic given by, 

m 

S m = y^(sj~ (s)) = -N oac (E m+ i) m = l,2,...n (3) 

i=l 

Once the unfolding is performed, on an average, a unit 
interval of the spectrum will have one level. Hence, S m 
represents the cumulative deviation until mth level, of 
the ith unfolded level from i. This quantity has a formal 
analogy with a time series. If the index i represents the 
(scaled) time, then (s, — (s)) represents the actual value 
assumed by the series at the ith time instant. Following 
Ref. H|, we take the power spectrum of S m as, 



S(f) = \6ff< 



1 



t 5> 



exp 



— 2mfm 

71 + 1 



(4) 



where / = 1, 2, . . . , n; Sf is the Fourier transform of S m . 
Now, we will present results to infer that (S(f)) oc /~ 7 
and 7 depends on the degree of chaos in the system. 

Before we plunge into the results, we show the statistic 
S m for our models. For the quartic oscillator, using the 
results in Refs. [l3j we symmetry decompose the level 
density to obtain the asymptotic integrated level density 
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FIG. 1: (Color online) S m — m curve for (a) kicked top and 
(b) quartic oscillator. The solid curve corresponds to chaotic 
case (k = 7,a = 20) and the dotted curve to the intermediate 
region (k — 3, a = 12), and the dashed curve to regular limit 
(k = l,a = 0). 



for the A\ representation, 
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where F(.) is the Gauss' Hypergeometric function. We 
use this expression to unfold the quartic oscillator levels. 
Fig. displays S m — m for the quartic oscillator and the 
kicked top for a choice of 3 parameters in the regular, 
chaotic and the transition region. In Fig. ^a) for the 
kicked top, the dashed curve corresponding to nearly the 
Poisson spectrum (k = 1) differs markedly from the solid 
curve for almost the GOE (k = 7) limit. We also plot an 
intermediate case (dotted curve) to show the transition 
taking place from Poisson to GOE type spectrum. This 
intermediate case at k = 3 is also qualitatively different. 
Each series displays slightly different memory effects cor- 
responding to various shades of anti-persistent time se- 
ries. In fact, such time series are known to display J -7 
noise and we expect similar result based on this analogy. 
We observe similar features for the quartic oscillator in 
Fig. njb) for a = 6,11.5,19.5, as reported in an earlier 
work of Bohigas • 

In Figs. [21 and |3 we display the ensemble averaged 
power spectrum of S m . Ensemble averaging is done as 
follows : for each a of the quartic oscillator we obtain 
2000 levels. After leaving out the first 200 levels, we 
create 3 sequences of 600 levels each. We further obtain 
similar sequences from more values of a separated by 
da = 0.1. For instance, a = 11.5 in Fig. CJb) corresponds 
to ensemble average in the range a = 11.2 — 11. 8 in steps 
of 0.1. Hence, the results (except for the integrable case 
in Fig. EUa)) represent an average over 21 level sequences 
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FIG. 2: Power spectrum of S n for kicked top at (a) k = 1, 
(b)fc = 2.0, (c) fc = 3.0 and (d) k = 7.0. The solid line is the 
least squares fit. The slope 7 is indicated in each graph. 

of length 600 each. For the integrable case, the ensemble 
consists 9 level sequences (3 from each of a = 0, 2 and 
6) of length 600 each. Similar averaging is done for the 
kicked top with 11 level sequences of length 800 each. 
For the regular and chaotic limits shown in Figs. |2Ia,d) 
and|3fa,d), there is a good agreement with the predicted 
slopes of 7 = 2, 1 respectively 0. There are deviations 
for / < f m [ n due to the effect of short periodic orbits 
with scaled time period / m i n . The deviation for large 
/ arises, partly, from approaching the Nyquist frequency 
at / = (71 + l)/2. It is clear from Figs. I2l3r b.c) that as 
the oscillator and the top explore the intermediate region 
between the regular and the chaotic limits, the slope 7 
smoothly changes from 2 to 1 . Intuitively, we can expect 
this because other statistic' in the transition regime, e.g, 
spacing distribution, generally vary smoothly too in this 
regime. 

In order to obtain a global picture, we compute the 
exponent 7 for a range of parameter values of a and k 
in the order-chaos transition region. In RMT, this tran- 
sition is described by a semi-empirical spacing distribu- 
tion P(s; (3) characterised by the parameter f3. As (3 is 
varied from to 1, P(s; (3) changes smoothly from Pois- 
son to GOE, reflecting the change in classical dynamics 
from regularity to chaos. Here, we use the distribution 
Pt(s]/3) due to Izrailev 0. For f3 = 0, Pt(s;/3) gives a 
Poissonian form (integrable limit) and for (3 = 1 it closely 
approximates the GOE distribution (chaotic limit). The 
intermediate values, < (3 < 1, correspond to order- 
chaos transition. Fig. 01a,b) shows that the series of (3 
and 7 display similar trends and are strongly correlated 
for both the top and the oscillator. It is known that the 
fraction of regular regions in phase space is correlated 
with a parameter like (3 that characterizes the change in 



FIG. 3: Power spectrum of S m for quartic oscillator at (a) 
a=0, 2,6 (b) a=7.5 (c) a=11.5 and (d) a=19.5 . The solid 
lines are the least squares fit with intercept shifted for clarity. 
The slope 7 is indicated in each graph. The arrows in (b-d) 
indicate log/ mm , the time period of the short periodic orbit. 



spacing distribution |16(. Hence we infer that the ex- 
ponent 7 in the power spectrum reflects the qualitative 
trends in the classical dynamics of the system. In gen- 
eral, we have shown that (S(f)} oc /~ 7 , where the value 
of 1 < 7 < 2 relates to the nature of classical dynamics 
for the oscillator and the top. 

As opposed to purely random matrices, the semiclas- 
sical systems deviate from /~ 7 scaling Q for / < / mm , 
where / mm is the period i mm of the shortest periodic 
orbit (PO) scaled by the Heisenberg time (tn — 2irh) [tJ, 
i.e, / mm = (n + l)imin/^ff- This corresponds to short 
POs being system specific features and in the correspond- 
ing large energy scales universality breaks down leading 
to deviations from RMT based results • Since the spec- 
tral form factor is linear only for times t/tjj < 1, which 
is necessary to realise 1// noise the actual range of 
scaling is restricted to / mm <C / <C fn- Since all times 
are scaled by tn, at t = tn, we have fH = {n + 1). 
The short PO in the quartic oscillator is the 'channel or- 
bit', as it is referred to in the literature, with the initial 
condition (qi — 0, Q2,Pi — 0,p2)- This can be identi- 
fied by taking discrete Fourier transform of the scaled 
energies of the oscillator ■ The period of this orbit is 
t(E) = E^ 1 / 4 ^Y{s/i) ■ Among all the sequence of levels 
that form the ensemble let the largest level be E = E'. 
Then, the period of the short PO is t m %n — t{E'). For 
instance, at a = 11.5, we have E' = 2579.05 (h = 1) 
and n — 600 and this provides the theoretical bounds 
for scaling to be log(/ TO j n ) = 1.54 and log(/#) = 2.77, 
as indicated in Fig. EUc). It is evident from Fig. Blb-d) 
that (S(f)) oc / 1 in almost the entire theoretically ex- 
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FIG. 4: The parameter f3 (squares) in Pr(s) as a function 
of chaos parameter for (a) the kicked top and (b) quartic 
oscillator. Also plotted is the shifted slope 2 — 7 (circles) of 
the power spectrum (S(f)). The slope is shifted such that 
y— axis lies in the range to 1. 



pected range f m \ n <C / <C fjj- This scaling range can be 
increased by probing deep in the semiclassical (h — > 0) 
regime, since tn ^ ft -1 ^ i? 3 ^ 4 for the quartic oscillator. 
But this leads to large eigenvalue problems that may not 
be computationally feasible at present. 

In Fig. \E\ we show the effect of short periodic orbits 
on the power spectrum of 5 m . Note that for the quartic 
oscillator at a = 20, 7 = 1.073 whereas for a — 30 we 
have 7 = 1.322 and (3 also shows a similar trend (see 
FigEJ- In general, we would expect that as a increases 
monotonically, chaos also increases and hence agreement 
with RMT should get better. But numbers quoted above 
show that, roughly speaking, chaos at a — 20 is more 
than at a = 30. This 'anomalous' feature is the effect 
of oscillating stability of the short PO. At a — 20,30 
the short PO undergoes an anti-pitchfork and a pitchfork 
bifurcation respectively accompanied by local changes in 
the phase space structure ^tJ- This, in turn, affects the 
spectral levels. It is known that the short POs influence a 
series of eigenstates, called the localized or sometimes the 
scarred states, in the quartic oscillator spectrum and they 
deviate strongly from RMT for eigenvector statistics [llj, 
[Lsl . If we remove the spacings that involve localized 
states, then we might expect the resulting distribution to 
show a better agreement with RMT. This is like removing 
the effect of short POs in the spectrum. The dotted 
line in Fig. 0is the usual power spectrum and the solid 
line is the one whose spacings involving localised states 
are removed. In this example the ensemble has just 3 
sequences of 600 levels each which is reflected in large 
amplitude of fluctuations. The power spectrum changes 
character for / < / m i n and the range of validity of power 
law gets better. This is a manifestation of the effect of 
short POs in the spectrum. 



FIG. 5: (Color online) The effect of short periodic orbits in the 
spectrum. (S(f)) for quartic oscillator at a — 30 including the 
localized states (dotted). The arrow indicates / mm . The solid 
line is the power spectrum after localized states are removed. 

In summary, we have shown that the spectral fluctu- 
ations in the quartic oscillator and the kicked top dis- 
play /~ 7 noise, where the value of 7 within the limits 
1 < 7 < 2 reflects the underlying nature of classical dy- 
namics, namely, regular or chaotic or a mixture thereof. 
We show the effect of short POs on the spectral fluctua- 
tions. We expect / _7 type noise in the level fluctuations 
to be an inherent characteristic of quantum systems. 

After this work was completed, qualitatively similar re- 
sults obtained on a billiard system were posted on LANL 
archive |2JJ|. We thank Prof. V. B. Sheorey and Dr. Dilip 
Angom for useful discussions. 
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